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ABSTRACT 

This  paper  presents  a  new  concept  for  the  remote  detection  of  radioactive  materials.  The 
concept  is  based  on  the  detection  of  electromagnetic  signatures  in  the  vicinity  of  radioactive 
material  and  can  enable  stand-off  detection  at  distances  greater  than  100  m.  Radioactive 
materials  emit  gamma  rays  which  ionize  the  surrounding  air.  The  ionized  electrons  rapidly 

attach  to  oxygen  molecules  fonning  07  ions.  The  density  of  07  around  radioactive  material 
can  be  several  orders  of  magnitude  greater  than  background  levels.  The  elevated  population  of 
07  extends  several  meters  around  the  radioactive  material.  Electrons  are  easily  photo-detached 
from  O7  ions  by  laser  radiation.  The  photo-detached  electrons,  in  the  presence  of  laser 
radiation,  initiate  avalanche  ionization  which  results  in  a  rapid  increase  in  electron  density.  The 
rise  in  electron  density  induces  a  frequency  modulation  on  a  probe  beam  which  becomes  a  direct 
spectral  signature  for  the  presence  of  radioactive  material. 

I.  INTRODUCTION 

Radioactive  materials  emit  gamma  rays  which  ionize  the  surrounding  air  producing  free 
electrons  which  rapidly  attach  to  oxygen  molecules  forming  07  ions.  The  density  of  negative 
ions  is  far  greater  than  the  free  electron  density  [1,2].  This  is  true  even  for  ambient  levels  of 
radioactive  material. 

A  recently  proposed  radioactivity  detection  concept  is  based  on  a  high  power  THz  pulse 
inducing  avalanche  breakdown  and  spark  fonnation  in  the  vicinity  of  the  radioactive  material  [3- 
6].  In  that  concept  a  THz  pulse  is  focused  to  a  volume  Vfocal  near  the  radioactive  material.  In  the 

absence  of  the  radioactive  material  the  electron  density  Ne  is  sufficiently  small  so  that 
Ne  Vfocal  «  1  .  In  this  case  the  probability  of  an  electron  appearing  in  the  focal  volume  during 
the  THz  pulse  is  negligible  and  the  breakdown  probability  is  very  low.  In  the  presence  of 
radioactivity,  however,  Ne  Vfocal  >  1  ,  which  will  lead  to  avalanche  ionization  in  the  volume 
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Vfocal  provided  the  THz  pulse  intensity  is  above  a  threshold  level  and  the  pulse  length  is 

sufficiently  long.  Breakdown  and  spark  formation  require  that  the  collisional  ionization  rate 
exceed  the  attachment  rate.  A  THz-based  radioactivity  detection  concept  may  have  unique 
advantages  depending  on  the  stand-off  distance  and  atmospheric  conditions. 

In  this  paper  we  propose  and  analyze  an  alternative  detection  concept  using  a  high 
intensity  laser  beam  and  probe  beam  to  detect  electromagnetic  signatures  in  the  vicinity  of 
radioactive  material.  Propagation  of  high  power  short  pulse  lasers  in  the  atmosphere  over 
extended  distances  has  been  analyzed  [7]  and  experimentally  characterized  [8-10].  Since  the 
negative  ions  produced  by  the  radioactive  material  have  a  low  ionization  potential  (0.46  eV)  they 
can  be  photo-detached  by  laser  radiation  (~  0.8  -lpm).  Our  detection  approach  is  based  on  a 
probe  radiation  beam  undergoing  a  frequency  modulation  while  propagating  in  a  temporally 
increasing  electron  density.  The  frequency  modulation  on  the  probe  beam  becomes  a  spectral 
signature  for  the  presence  of  radioactive  material.  A  schematic  of  the  detection  concept  is  shown 
in  Fig.  1. 
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Figure  1 .  Schematic  of  active  remote  radioactivity  detection  concept.  Laser  radiation  photo- 
detaches  electrons  from  Of  ions,  providing  electrons  for  an  avalanche  (collisional)  ionization 
process  which  increases  the  electron  density.  A  probe  beam  in  the  presence  of  a  temporally 
increasing  electron  density  undergoes  a  frequency  modulation  which  becomes  a  spectral 
signature  of  radioactivity. 

II  RADIOACTIVITY  DETECTION  CONCEPT 

In  the  following  subsections  the  various  elements  of  the  detection  concept  are  discussed 
and  analyzed.  Details  of  the  analytical  model  are  presented  in  the  appendix. 

a)  Radiation  Enhancement  Factor 

A  gamma  ray  propagating  through  matter  can  interact  through  several  processes, 
including  Rayleigh  scattering,  photoelectric  and  Compton  effects,  pair  production,  etc  [1 1].  In 
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air,  photoelectric  absorption  dominates  at  low  photon  energies  (<  25  keV)  while  at  high  energies 
(~  25  keV  -  3MeV)  Compton  processes  dominate.  As  the  gamma  ray  propagates  in  air  it  loses 
energy  in  a  cascading  process  and  its  mean  free  path  L  decreases.  Gamma  rays  having  energies 

of  IMeV,  50keV,  25keV  or  lOkeV,  have  mean  free  paths  in  air  of  Ly  «  130m ,  38m,  19m  or 

1 .9m,  respectively.  Note  that  values  for  the  range  indicated  here  apply  to  the  case  where  a 
gamma  ray  with  the  given  energy  is  completely  absorbed  in  a  single  interaction  (collision).  The 
actual  mean  free  path  for  a  given  initial  energy  is  somewhat  less  than  that  for  a  single  interaction. 

The  ionization  rate  due  to  background  (ambient)  radioactivity  is  {cl  Ne  /  dt)amh  =  Qrad .  At 
or  near  ground  level,  the  background  ionization  rate  is  typically  in  the  range 
Qrad  ~  1 0  —  30  pairs/(cm3  -  sec)  .  The  gamma  rays  emitted  by  radioactive  material  ionize  the 
surrounding  air.  In  the  presence  of  radioactive  material  the  ionization  rate  (due  to  only  radiation) 
can  be  greatly  enhanced  by  a  factor  arad  »  1  and  (d  Ne  /  dt)rad  =  amd  Qmd  .  Various  air 

chemistry  processes  such  as  attachment  and  recombination  limit  the  electron  density  and  are 
discussed  later  and  in  the  appendix. 

The  enhancement  factor  amd  can  be  estimated  as  follows.  For  a  radioactive  material  of 
mass  Mmd  the  number  of  disintegrations  per  second  is  vmd  =  Mmd  Amd ,  where  Arad  is  the 

specific  activity  associated  with  the  material.  For  example,  for  60 Co  , 

Amd  =  l.lxlCfCi/g  =  4.1xl0lj  disintegration/^- sec)  ,  where  a  Ci  (Curie)  is  3.7xl010 

disintegrations/sec.  In  the  case  of  60Co  each  disintegration  results  in  two  gammas  of  energy 
Ev  =  1.173MeV  and  Ev  =  1.332MeV  which  have  a  range  in  air  of  ~  130m .  In  air  the 
high  energy  gammas  generate  high  energy  electrons,  via  Compton  and  photoelectric  processes, 
which  undergo  a  cascading  process  to  sufficiently  low  energy  to  attach  to  02  molecules  forming 

0~  ions.  In  the  cascading  process  the  electrons  lose  an  amount  of  energy  A E  «  34  eV  per 
collision  in  air  which  results  in  both  ionization  and  electronic  excitation  (~  equal  amounts  of 
energy  lost  to  each  process)  [12].  A  high  energy  electron  with  energy  Ee  therefore  generates 

~  Eel  A E  low  energy  electrons.  In  the  case  of  Compton  absorption,  the  maximum  electron 

energy  is  Ee  =  (2ay  /(\  +  2ay))  C,,  max  where  ay  =  Eymax  /me2 .  For  example,  a  1  MeV  gamma 

ray  in  air  generates  Compton  electrons  having  a  maximum  energy  of  Ee  =  0.8  MeV  and  average 

energy  of  {EeJ  =  0.44  MeV .  It  should  be  noted  that  the  range  of  high  energy  electrons  is  much 

less  than  the  range  of  the  high  energy  gammas.  For  example,  an  electron  having  an  energy  of  1 
MeV,  0.5  MeV  or  0.1  MeV  has  a  range  in  air  of  4.6  m,  3  m,  or  2.1  m,  respectively.  Here,  again 
as  with  the  gamma  rays  mentioned  earlier,  the  above  electron  ranges  corresponds  to  the  case 
where  an  electron  is  completely  stopped  in  one  interaction  (collision). 

For  a  small  spherical  source  of  radioactivity  the  steady  state  density  of  emitted  gamma  rays  is 


N„  =  exp(-  R/ L  ) , 


(1) 


4  ncR~ 

where  R  is  the  distance  from  the  radioactive  material,  L  ,  is  the  effective  range  (mean  free  path) 
of  the  gamma  rays  in  air  which  is  a  function  of  the  gamma  ray  energy,  E  ,  and  k  is  the  number 
of  gammas  emitted  per  disintegration.  The  rate  of  change  of  electron  density  is  given  by 
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where 


— ~  ~  (amd  + 1 ) Qmd  +  a'r  chemistry  and  ionization  terms, 
dt 


(2) 


a 


rad 


»  C(<J 


)NairNy 


A  E  Qrad 


VradKy  (ge)  1  CXP(~ /?  / A,  ) 

4^(V,}  AE  Qrad  R2 


(3) 


is  the  radiation  enhancement  factor.  An  equivalent  result  defining  the  detectable  mass  of  the 
radioactivity  material  as  a  function  of  the  various  parameters  was  obtained  in  Ref.  6. 

In  the  absence  of  radioactive  material  amd  =  0 .  In  Eq.(3)  (<?r-e)  =  (crr-e)CJ  +  (°V-e)P£, 
is  the  effective  average  cross  section  for  electron  generation  by  gammas,  i.e.,  Compton 
absorption  and  photoelectric  processes,  (C7_e)  =  ((crr_e)  AC.)  is  the  average  mean  free  path  for 

electron  generation  by  gammas  and  Nair  =  2.7  x  1019  cm  is  the  air  density  (Loschmidt’s  number) 
at  STP.  The  cross  section  <j r_e(E)  is  a  function  of  the  gamma  ray  energy  and  the  average  cross- 

section  is  given  by  =  J  dE  f  (E)ay_JE) ,  where  / (E)  is  the  gamma  ray  energy 

distribution  function.  Over  the  range  of  energies  ~  0.05  -  IMeV  Compton  absorption  dominates 
and  (<Jy_}]cA  ~  1.4xl0-24  cm2 ,  while  for  energies  below  ~  0.05  MeV  the  photoelectric  effect 

dominates  with  ^ crr_e^pE  increasing  from  ~1.4xl0-24  cm2  at  0.05  MeV  to  ~2xl0_22cm2  at  0.01 

MeV.  In  Fig.  2  the  radiation  enhancement  factor  amd  is  plotted  as  a  function  of  the  distance 

from  the  radioactive  source,  R  .  This  plot  is  for  samples  containing  1  mg  and  lOmg  of  60 Co  and 
indicates  that  the  enhanced  level  is  significant  for  ranges  extending  up  to  several  meters.  As  an 
example,  for  (Ee)  =  0.5MeV,  Mmd  =  lOmg,  vrad  =  MradArad  =  8.2xlOu  disintegrations/sec , 

Ky  =  2  ,  R  =  50cm  ,  (Lr  }j  =  100m  and  Qmd  =  20  disintcgrations/(cnr!- sec)  the  radiation 

enhancement  factor  is  amd  w  2xl06  which  is  far  above  the  background  level. 


Figure  2.  Radiation  enhancement 
factor  plotted  as  a  function  of 
distance  from  the  radioactive  source 
for  1  mg  and  10  mg  of  60 Co  . 
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b)  Electron  and  Ion  Density  Evolution  (Air  Chemistry) 

To  determine  the  frequency  modulation  on  a  probe  pulse  it  is  necessary  to  follow  the  time 
evolution  of  the  electron  and  negative  ion  density,  which  are  sensitive  functions  of  air  chemistry 
processes  [13]  and  electron  heating  by  the  laser  radiation.  The  source  terms  for  the  electrons 
include  radioactivity,  detachment,  photo-detachment  and  photo-ionization,  while  the  loss  terms 
include  various  attachment  and  recombination  processes  including  aerosols.  The  expressions  for 
the  rate  of  change  of  electron  density  N e  and  negative  ion  density  N _  are  [14-17], 


dNJdt  =  (\  +  ccrad)Qrad  +  Se  -  Le, 

(4a) 

dNJdt  =  S_-  L_, 

(4b) 

where  Se  represents  the  various  electron  source  terms,  Le  is  the  electron  loss  terms,  S_ 
represents  the  ion  source,  L_  is  the  ion  loss  tenns.  The  various  source  and  loss  terms  are 
discussed  in  the  appendix. 

The  effect  of  radioactivity  is  represented  by  the  first  tenn  on  the  right  hand  side  of  the 
electron  rate  equation,  Eq.  (4a).  The  steady  state  electron  and  negative  ion  densities  are  given  by 

Ne  «  (/?„  Nn  /  ?/)  V(1  +  a  rad  )Qmd  /  P+  +  <\+<Xrad)Qrad  1  *1  *  iPn  Nn  /  Tj)J  (1  +  U  J)  Qrad  t  P+  ,  and 

+  a rad) Qrad  / P+  [4],  where  N„  is  the  neutral  air  density  ( Nn  ~  Nair  for  low  levels  of 
ionization),  P  »  2x10  6  cm3 /sec  is  the  recombination  rate,  //  ~  10s  sec  is  the  attachment 
rate  and  pn  ~  (5  - 10)  x  10  19  cm3  /  sec  is  the  negative  ion  detachment  rate  due  to  collisions  with 
neutrals  (see  appendix  for  details).  In  the  absence  of  radioactive  material  ( arcd  =  0  )  the  ambient 
(background)  electron  and  negative  ion  density  are  Ne  «  10  cm  and  N_  ~  3x10'  cm”3, 
where  Qrad  =  20  cm  sec  ”1 .  At  a  distance  of  4  m  from  a  radioactive  source  containing  10  mg  of 
60 Co  ,  the  radioactive  enhancement  factor  is  arad  =2.2xl04  and  the  elevated  electron  and 
negative  ion  densities  are  Ne »  0.2  cm-3  and  N_  »  4.7 xlO5  cm”3 . 

The  various  source  and  loss  terms  in  Eqs.  (4a, b),  in  particular  the  collisional  ionization 
rate,  are  functions  of  the  electron  temperature  (see  appendix).  The  electron  temperature  is 
determined  by  the  collisional  electron  heating  (Ohmic  heating)  by  the  laser  radiation  and  the 
cooling  effect  resulting  from  excitation  of  vibrational  modes  of  the  air  molecules.  The  equations 
for  the  electron  temperature  Te  is  [18], 

(3 / 2)  8 ( W,  r, ) / 0 <  =  ( J •  E)  +  (3 / 2) ( W,  / r,^ ) (r,  - r„ )  -  (/„  0  AT,  / 0 < ,  (5) 

where  Te  is  the  electron  temperature,  (j  ■  E)  is  the  Ohmic  heating  rate,  zcool  is  the  electron 
cooling  time  due  to  inelastic  collisions  (see  appendix),  Teo  =  0.025  eV  is  the  ambient  electron 
room  temperature,  Nair  is  the  background  air  density  and  Uion  is  the  effective  ionization 
potential  of  air  (  ~  10  eV).  The  heating  and  cooling  terms  on  the  right  hand  sides  of  Eq.  (5)  are 
discussed  in  the  appendix. 

c)  Frequency  Modulation  on  a  Probe  Beam 
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A  probe  beam  propagating  through  a  region  of  space  in  which  the  electron  density  is 
changing  with  time  will  undergo  a  frequency  change.  The  electron  density  in  the  vicinity  of  the 
radioactive  source  and  under  the  influence  of  the  laser  radiation  varies  in  space  and  in  time. 
Consequently  the  frequency/wavenumber  of  an  electromagnetic  probe  beam  propagating  in  such 
a  plasma  will  vary  in  space  and  in  time.  The  one-dimensional  wave  equation 
(02  / dz2  -  c~2d2  / dt2)A(z,t )  =  c~2a>2(z,t)A(z,t )  can  be  used  to  determine  the 

frequency/wavenumber  shift.  Here,  co  (z,t)  =  (4 tt q2Ne(z,t)/ in)'  2  is  the  plasma  frequency  and 
A(z,  t )  is  the  vector  potential  associated  with  the  probe.  The  vector  potential  can  be  expressed  in 
tenns  of  slowly-varying  amplitude  and  phase,  in  the 

form  A(z,t)  =  (l/2)B(z,t)exp[i(k0z-a>0t)  +  i  9(z,t)]  +  c.c.,  where  coo  is  the  incident  probe 
frequency  and  9{z,t )  is  the  phase.  The  frequency  of  the  probe  is  0Jprobe(z,t)  =  a>0  -  d0(z,t)/dt , 
where  the  phase  6  satisfies  (d/dz  +  c~'d/ dt)0(z,t)  =  -co2(z,t)/(2a>0c) .  The  solution  is 

rj<L 

=  -(2coa  c)~l  jd//co2(//,r  +  i//c)  where  rj  =  z,  r  =  t  —  zlc  and  L  is  the  probe  beam’s 

o 

interaction  length.  The  frequency  shift  on  the  probe  is  A co(tj,  t)  =  —  d0(rj, r)l dr  .  In  particular 
if  the  electron  density  is  increasing  with  time  the  probe  frequency  will  be  increased.  For  a 
spatially  uniform,  temporally  varying  plasma  density,  the  frequency  shift  on  the  probe  beam  is 
A co(z,t)  =  (2cooy\(D2p(t)  -  co2p(t-z/c)),  (6) 

where  0  <z<L. 

As  an  illustration  we  consider  the  case  where  the  rise  in  electron  density  is  exponential  in 
time  and  spatially  uniform  within  a  region  0  <  z  <  L  .  Here  the  plasma  frequency  is  given  by 

co2(z,t )  =  orpo  exp (vion  t)g(z)  ,  where  vion  is  the  effective  ionization  rate  and  g(z)  =  1,  0  <  z  <L  , 
defines  the  axial  extent  of  the  plasma  region.  The  frequency  modulation  on  the  probe  beam  is 
given  by  coprobe(z,t)  =  (Oo  +  Aa>(z,t) ,  where 

A a>(z,t)  =  (o2po  /(2coJcxp(viim  t)(l  -  exp(-i/f)n  z/c)).  (7) 

The  maximum  fractional  frequency  shift  occurs  for  z  >  c  /  vion  ~  1  cm  and  is 

A®max  /  (0„  =  (OJl0  /  2<xr)  exp {vion  t ) .  The  effective  ionization  rate  can  vary  widely  but  is  typically 

v ion  ~  1011  sec  1 .  The  frequency  modulation  on  the  probe  beam  can  be  substantial  and  is  readily 

measurable. 

III.  RADIOACTIVITY  DETECTION  EXAMPLE 

To  illustrate  the  radioactivity  detection  concept,  i.e.,  the  frequency  modulation  induced 
on  a  probe  beam  by  ionizing  the  air  near  the  radioactive  material,  we  consider  the  following 
example.  In  the  numerical  examples  the  system  quantities,  i.e.,  electron  density,  ion  density, 
electron  temperature  and  frequency  modulation  on  a  probe  laser,  are  calculated  both  in  the 
absence  ( amd  =  0  )  and  in  the  presence  of  radioactive  material  by  numerically  solving  Eqs  (4-6). 
The  radioactive  source  is  assumed  strong  enough  to  produce  a  radiation  enhancement  factor  of 
amd  =  103 .  The  radiation  enhancement  factor  value  is  consistent  with  low  quantities  of 
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Electron  Density  N<\  l|  c|rf|  Electron  Density  \\  l|  cfrf| 


radioactive  material  as  shown  in  Fig.  2.  We  take  the  ionizing  laser  to  have  a  peak  intensity  of 
I  k  -  160  GW/cm2  and  pulse  duration  of  rlaser  =  lnsec  .  In  these  examples,  the  probe  beam  is 
taken  to  be  a  millimeter  wave  source  of  frequency  fprobe  =  94  GHz  ,  (A  be  =  3.2  mm).  The 
critical  electron  density,  associated  with  the  probe  frequency, 

M probe  =  (0p.cn,  =  5.64  x  1 04  «''f([cin  3],  is  necrit  =  1014  cm  3 .  The  background  radiation  is  taken 
to  be  Qmd  =  30  disintegrations /(cm’  -sec) 


Figure  3  a  Electron  density  as  a 
function  of  time  in  the  absence  of 
any  external  radioactivity 
( amd  =  0  )•  The  laser  parameters 
are, 

1  =  lpm,  I peak  =  160  GW/cm2 

and  Z  laser  =  lnSCC' 

1  2  3  4  5 

Time  |  njancQec| 


Time  |  njanoQec| 


Figure  3b.  Electron  density  versus 
time  in  the  presence  of 
radioactivity.  The  laser  parameters 
are  the  same  as  in  Fig.3a. 


In  the  absence  of  radioactivity  the  ionizing  laser  intensity  is  just  below  the  breakdown  level,  i.e., 
the  electron  density  is  low,  and  there  is  virtually  no  frequency  modulation  on  the  probe  beam. 
Figure  3b  shows  the  electron  density  as  a  function  of  time  in  the  presence  of  radioactive  material 
( amd  =  103).  The  electron  density  at  the  end  of  the  ionizing  laser  pulse  approaches  the  value  of 

ne  =  101 3  cm  which  is  an  order  of  magnitude  less  than  the  critical  electron  density. 
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The  frequency  modulation  on  the  probe  millimeter  wave  beam  is  shown  in  Fig.  4.  In  the 
absence  of  radioactive  material  there  is  no  frequency  modulation  on  the  probe.  However,  for 
amd  =  103  the  fractional  frequency  modulation  is  significant  and  equal  to  ~  5  % ,  which  is 
readily  detectable. 


Figure  4.  Fractional  frequency  shift 
A  col  0)o  [%]  versus  time  in  the 
presence  of  radioactive  material 
°md  =  103  at  the  probe  interaction 
distance  of  L  =  10cm.  The  laser 
parameters  are  the  same  as  in  Fig.  3a. 


The  fractional  frequency  shift  on  the  probe  as  a  function  of  both  axial  interaction  distance  L  and 
time  is  shown  in  Fig.  5. 
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Figure  5.  Fractional  frequency  shift 
Aoj/  (oo  [%]  versus  time  and  probe 
interaction  distance  L  in  the 
presence  of  radioactive  material 
( amd  =  103).  The  laser  parameters 
are  the  same  as  in  Fig.  3  a. 
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IV.  DISCUSSION  AND  SUMMARY 

In  this  paper  a  concept  is  proposed  and  analyzed  for  active  remote  detection  of 
radioactive  materials.  The  enhanced  levels  of  ionization  associated  with  the  presence  of 
radioactivity  may  be  detected  using  a  combination  of  an  ionizing  laser  and  probe  beam.  The 
detection  concept  depends  on  the  fact  that  the  free  electrons  produced  by  gamma  rays  emitted  by 
the  radioactive  material  attach  rapidly  to  oxygen  molecules  fonning  negative  ions  0~ .  In  the 
vicinity  of  radioactive  material  the  density  of  negative  ions  is  elevated  compared  to  the 
background  density,  while  the  free  electron  density  remains  extremely  low.  The  negative 
oxygen  ions  can  be  readily  photo-detached  with  laser  radiation.  The  photo-detached  electrons,  in 
the  presence  of  the  laser  radiation,  initiate  an  avalanche  ionization  process  which  greatly 
increases  the  electron  density.  To  detect  this  change  in  electron  density  a  probe  beam  is 
employed.  The  probe  beam  undergoes  frequency  modulation  in  the  presence  of  the  temporally 
changing  electron  density.  The  frequency  modulation  on  the  probe  beam  is  a  signature  of  the 
presence  of  radioactive  material.  Our  analysis  indicates  that  a  measurable  frequency  shift  can  be 
expected  for  relatively  small  amounts  of  radioactive  material. 

Plans  are  underway  to  experimentally  demonstrate  this  detection  concept  at  the 
University  of  Maryland.  The  experiments  employ  a  commercial  source  of  alpha  particles  to 
simulate  the  nuclear  material.  Note  that  in  an  actual  detection  device  gammas  not  alphas  would 
provide  the  initial  negative  ions.  The  nuclear  material  is  210 Po ,  producing  5  MeV  alpha 
particles.  Using  a  radioactive  source  of  alphas  in  place  of  a  gamma  ray  source  avoids  the  safety 
issues  associated  with  radioactive  material.  Note  that  5  MeV  alphas  have  a  short  range  (~ 

3.5  cm  )  in  air.  The  commercial  source  produces  20  mCi  of  radioactivity  from  210  Po  inside  an 
open  metallic  tube.  Each  5  MeV  alpha  will  generate  ~  1.4  x  105  electrons  which  attach  to  02 
fonning  02  ions.  The  ()2  ions  are  blown  out  of  the  tube  and  are  photo-detached  by  laser 
radiation. 

The  photo-detaching  and  ionizing  beam  in  the  experiments  is  a  Nd:YAG,  5GW,  100  psec 
laser.  The  diagnostics  include  visible  and  extreme  ultraviolet  spectrometers,  interferometers  for 
measuring  electron  density  profdes,  microwave  and  ultrafast  pulse  diagnostics.  The  probe  beam 
used  in  these  experiments  is  a  30  GHz,  circularly  polarized  continuous  wave  microwave  source. 
Focusing  the  laser  radiation  near  the  210 Po  source  photo-detaches  the  negative  ions  and  initiates 
an  avalanche  ionization  process.  The  temporal  rise  in  electron  density  produces  a  detectable 
frequency  modulation  on  the  probe  beam. 

The  purpose  of  this  paper  is  to  present  a  new  concept  for  the  detection  of  radioactive 
material.  Our  analysis  indicates  that  a  measurable  frequency  modulation  on  a  probe  beam  is 
induced  even  in  the  presence  of  small  quantities  of  radioactive  material.  A  more  rigorous 
analysis  requires  a  more  detailed  model  of  the  air  chemistry  processes. 
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APPENDIX 


Electron  and  Ion  Density  Dynamics 

The  rate  of  change  of  electron  density  Ne  and  the  negative  ion  density  N  ,  Eqs.  (4a)  and 


(4b),  are  given  by 

dNe/dt  =  aradQrad  +  Se-Le,  (Al) 

dNJdt  =  S_  -  L_,  (A2) 

and  the  source  and  loss  terms  are  given  by 

Se  =  vcoll  Ne  +  vopt  N_  +  N_  Ne  +  pn  Nn  N_ ,  (A3) 

4  =  Pe  N+  Ne  +>lNe+  pAe  N A  Ne ,  (A4) 

S.  =  rjNe,  (A5) 

L  =  P+  N+N_  +  pn  NnN_  +  voptN_  +  P_N_Ne  +  pA_NAN_,  (A6) 


and  the  positive  ion  density  is  given  by  N+  =  Ne  +  N_  (charge  neutrality).  In  air,  the  negative 
ion  population  is  mainly  Oj .  In  these  equations  Nn  is  the  neutral  air  density,  N  4  is  the 
density  of  aerosols,  vcotl  is  the  collisional  ionization  rate,  v  t  is  the  photo-detachment  rate, 

Pe  ~  3xl0_x  cm'  sec-1  is  the  electron-ion  dissociative  recombination  rate, 

Pn  =  (5-10)xl0“19  cm3  sec-1  is  the  negative  ion  detachment  rate  due  to  collisions  with  neutrals, 
r]  ~  108  sec”1  is  the  electron  attachment  rate,  /?+  ~  2xl0-6  cm3  sec-1  is  the  negative  ion 
recombination  rate,  PA_  is  the  aerosol-ion  attachment  rate  and  pAe  is  the  aerosol-electron 
attachment  rate,  vcoll  is  the  collisional  ionization  rate,  and  P  is  the  negative  ion  detachment  rate 
due  to  collisions  with  electrons.  The  various  attachment  and  recombination  rates  are 
complicated  functions  of  the  system  parameters,  in  particular  the  electron  temperature.  In  our 
model  we  approximate  the  variouse  attachment  and  recombination  rates  as  constants  given  by 
the  above  values. 

Photo-Detachment  of  0~ 

The  ionization  potential  of  0~  is  0.46  eV  and  therefore  can  undergo  single-photon  photo¬ 
detachment  with  laser  radiation  of  wavelength  X  =  lpm  (1.24eV)  or  X  =  0.8pm  (1.55 eV).  The 
photo-detachment  rate  is  vopt  =  cropt  c  N  .  =  <Jopt  I0lhco,  where  c  N  h  -  I  o!  fico  is  the  incident 
photon  flux,  Io  is  the  laser  intensity  and  aopl  is  the  photo-detachment  cross  section.  The 
experimental  value  for  the  single-photon  ionization  cross  section  of  Cf"  is 
(jopt(X  =  \\im)  «  4.5xl0”19cnr  and  cropt (T  =  0.8 pm)  «  7.5x1 0”19  cm2  [16].  The  single-photon 
ionization  rate  for  0A  is  therefore, 

,  ,  [2.3,  2  =  1  pm 

v,[^]-/.[WW]|3j  (A7) 

For  the  parameters  considered  here  multi  photon  ionization  of  air  can  be  neglected  since  the 
ionization  potential  for  clean  air  ( 02  )  is  —  12  eV  and  requires  -8-10  photons. 
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Collisional  (Avalanche)  Ionization  of  Air 

oo 

The  collisional  ionization  rate  is  given  by  vcoll  =  N air  J  cl  v  <xco//(v)  v  /e(v) ,  where  acoll  is 

v7 

the  collisional  ionization  cross  section  and  fe(y)  is  the  electron  distribution  function  is  taken  to  be 
a  Maxwellian,  with  J  J  v  fe=  1 .  The  cross  section  can  be  approximated  near  the  threshold  by 
<jcoll  «  C(£-Uj)  where  s  =  m  v2/ 2  is  the  electron  energy,  U1  is  the  ionization  energy, 
v7  =  (2 Uj  / /u)12  and  C  is  a  constant  associated  with  the  molecule.  The  collisional  ionization  rate 
in  air  can  be  written  in  the  form  vcoll  =  vcoU(N2)  +  vcoll(02)  where  the  ionization  rates  of  N2  and 
02  have  the  form 

VcoiiiX)  =  vx(TJUxy2{UxITe  +  2)exp  (~UX/Te),  (AS) 

with  UNi  =15.6  eV,  C/0i  =12.1  eV  and,  at  STP,  vNi  =  7.6x10“  sec-1  and  vQi  =  10“  sec-1  [18]. 

Electron  Cooling  Rate 

The  various  inelastic  electron  cooling  processes  in  air  are  complicated  and  include 
excitation  of  vibrational,  rotational  and  electronic  states.  At  low  electron  energies  the  dominant 
inelastic  cooling  processes  is  the  excitation  of  vibrational  states  of  N, .  As  the  electron 
temperature  increases  the  electron  cooling  rate  decreases  due  to  collisional  excitation  of  the 
nitrogen  molecule  vibrational  states.  We  model  the  inelastic  electron  cooling  rate  as 
1/  rcoo/  =  (2/3)(2m/M eff)ve ,  where  Meff  is  the  effective  molecular  mass  [18]  and  ve  is  the 
electron  collision  frequency. 

Electron  Heating 

The  collisional  electron  heating  tenn  (J-E^>  in  Eq.  (5)  is  obtained  by  representing  the 
laser  field  and  electron  current  density  by  E  =  E0(z,t)exp(i X¥(z,t))ex  / 2  +  c.c.  and 
J  =  J0(z,t)exp(i'¥(z,t))ex  / 2  +  c.c.,  respectively,  where  'F  =  co(zlc-t) , 

J0(z,t )  =  i q2  Ne(z,t)(E0(z,t)/m)/(co  +  ive(z,t)) ,  EJV/cm]  =  27.4/]  2[W/cm2]  is  the  laser 
field  amplitude,  I0  =  c  \Eo |  / 8 n  is  the  intensity,  <x>  is  the  laser  frequency,  v(,  is  the  electron 
collision  frequency,  ex  is  a  unit  polarization  vector  along  the  x  -  axis  and  c.c.  denotes  the 
complex  conjugate.  The  electron  current  density  and  field  are  related  by 
d  J/dt  +  J  =  (®2/4;r)E,  where  cop  =  (4  nq2  Ne  /  m)112  is  the  plasma  frequency.  The 
electron  collision  frequency  is  the  sum  of  the  electron-neutral  and  electron-ion  contributions, 

=  v,„  +  E, ,  where  ujsec  ']  »  10  7  Nn  [cm  5]7;,|/2[eV]  and 
u^.fsec  ']  ~  1 0  5  N:,  [cm  ']  f  ’  2[eV] .  The  time-average  rate  of  change  of  electron  energy 
density  is  (J  •  Ej  =  (co2p  l%7t)E2eff  / vi; ,  where  Eeff  =  (ve  /  co)K\  +  v2  /  of )'  2  Eo  is  the  effective  laser 
electric  field.  In  practical  units 

( J  •  E>  [eV  -  cm-3  /  sec]  =  1.9x10' “3  Ae[cm -3]A2[nm]/0[W/cm2]ve[sec' lV(\+v2e  I  of) .  (A9) 

Typically  for  lasers  oo  »  ve  while  for  millimeter  waves  oo  can  be  comparable  to  ve . 
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